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Parametrization of the feedback Hamiltonian realizing a pure steady state 
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Feedback control is expected to considerably protect quantum states against decoherence caused 
by interaction between the system and environment. Especially, Markovian feedback scheme devel- 
oped by Wiseman can modify the properties of decoherence and eventually recover the purity of 
the steadystate of the corresponding master equation. This paper provides a condition for which 
the modified master equation has a pure steady state. By applying this condition to a two-qubit 
system, we obtain a complete parametrization of the feedback Hamiltonian such that the steady 
state becomes a maximally entangled state. 
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I. INTRODUCTION 

The time evolution of a quantum state under Marko- 
vian decoherence is usually modeled by a master equation 
in the Lindblad form 0, 



dp 



M 



(1) 



k=l 



T>[b k }p := b kP b\ - \b\b k p-\pb\b 



where H = W denotes a driving Hamiltonian, and b k are 
the Lindblad operators that generate decoherence. The 
dynamics Q often reduces the purity of p, where the 
purity is denned by p := Tr (p ). For example, a master 
equation that preserves the identity operator is strictly 
purity decreasing. (Conversely, any finite-dimensional 
purity-decreasing dynamics preserves the identity oper- 
ator 0|.) Since loss of purity is the most critical issue 
in many quantum information technologies Q, we must 
exploit effective methodologies to prepare desirable pure 
states even under decoherence. 

Quantum control shows a high possibility to overcome 
the problems of decoherence. In particular, quantum fil- 
tering theory, which was pioneered by Belavkin 0, 
enables us to feed measurement data continuously back 
in order to control quantum systems. Actually, in the 
framework of Markovian feedback proposed by Wiseman 
0, a data-dependent feedback Hamiltonian is continu- 
ously added to the system, and the structure of the cor- 
responding master equation is modified. Markovian feed- 
back method has been widely used for several types of 
quantum systems. For example, a master equation de- 
scribing a two-level atom was modified by the feedback 
so that any pure state becomes a steady state of the dy- 
namics 0. Further, in the case of a system consisting 
of two two-level atoms (qubits) coupled to a single-mode 
cavity field that is heavily damped, it was proved that 



the entanglement of the steady state was increased by 
the application of Markovian feedback [g ■ 

However, the selection of the data-dependent feedback 
Hamiltonian is usually based on intuitive observations. 
In other words, there are no systematic guidelines to 
modify the master equation to a desirable one. Actu- 
ally, the above-mentioned steady state of the controlled 
two-qubit system was not a maximally entangled state 
but a mixed state. 

Therefore, for finite-dimensional quantum systems, 
this paper characterizes a feedback Hamiltonian com- 
pletely such that the modified master equation has a 
pure steady state. By applying this condition to the two- 
qubit system, we obtain a parametrization of the feed- 
back Hamiltonian that causes the steady state to be not 
only a pure state but also a maximally entangled state. 



II. A CHARACTERIZATION OF THE 
FEEDBACK HAMILTONIAN 



A. The condition for a pure steady state 

The steady state of the master equation JTJ is a quan- 
tum state p satisfying Lp — 0. Remarkably, the steady 
state p is stable in the following sense. Let us de- 
fine the matrix elements pij (I < i,j < n) of an n- 
dimensional state p as follows. Set pa := vi, vi S K 
for 1 < i < n — 1, which leads to p nn = 1 — J2i v %- For 

For i < j 
An, and 



i > j set pij := Xij + ifMj with Ay, p^ € K. 
set pij — p*j . Collecting the real numbers v t , , Hj , 



Pij into a vector x £ 



we can describe the master 
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equation as an afiine equation x — Ax + a, where the 
matrix A and the vector a are uniquely determined by 
the above rule specifying x. Then, all the eigenvalues of 
A have non-positive real parts. This implies that for an 
initial state p + Sp {8p implies a small perturbation from 
p ), we can always find p(t) sufficiently close to p Q for all 
t > 0. Especially in the case when all the eigenvalues of 
A have negative real parts, p becomes a unique steady 
state corresponding to x Q = —A~ l a in the vector form; 
eventually every trajectory that obeys the dynamics must 



2 



converge into p . 

By combining the above discussion with the fact that 
pure states are the basis of many quantum information 
technologies, we aim to characterize a class of the Lind- 
blad operators bk and the Hamiltonian ff such that the 
master equation Q has a pure steady state. 

Theorem 1. The master equation 0) has a pure 
steady state if and only if bk (k = 1, . . . , M) and iff + 
(1/2) Sfc kfcfrfc have a common eigenvector \<f>). Then, the 
steady state is given by p = \4>){<j)\. 

Proof: Substituting p = \4>){<j)\ into the equation 
Cp = and multiplying by \cf>) on both sides, we have 



M 



= 0. 



(2) 



k=l 



Since the Schwarz inequality yields 

(#l& fc |MM) > \(4>\b k \4>)\\ 

Eq. @ is attained if and only if bk\(f>) — Pk\<t>) (k = 
1,...,M) holds for some p k € C. Thus, owing to 
b k \<j)){(j)\bl = \p k \ 2 \4>){<j>l the equation C\<f>){</>\ = is 
transformed to 



where 



F 



F|^)(^| + |0)^|Ft=O 

M , 



(3) 



k=l 



I MA I 



By denoting the /cth elements of |0) and F\<j>) by ^ and 
fk, respectively, we observe that the (j, k) element of the 
matrix equation is fj<fi% + <frjfk ~ 0' which becomes 



(k) + (it) 



fk 



0. 



when (j)j ^ and ^ 0. Then, /3j 



(4) 

fj/fa sat- 
isfies (3j + f3* — 0. Thus, /3j is a pure imaginary 
number: = i/3'., /3J S M. Therefore, Eq. gj) yields 
P'j — f3' k — 0. This leads to the conclusion that (3j is 
independent of the index j. Accordingly, we obtain 
fj = (3(f) j by defining a constant (3 := /3j, v j. When 
4>j = 0, it immediately yields fj = 0. Consequently, we 
obtain F\</>) — f3\cf>). It is easy to verify the "if condition. 

The proof of Theorem 1 is very similar to the one 
found in H: the theory of decoherence-free subspace 
(DFS) 111 ITU . We now consider a relationship 
between DFS and a pure steady state. Let the sys- 
tem Hilbert space TL decompose into a direct sum as 
TL = Hdf © Hj f , and partition the system state, the 
Hamiltonian, and the Lindblad operators into blocks as 
follows: 



P 



Pi P2 
P\ P3 



H 



ffl ff2 

Hi H a 



P k Qk 

Rk Sk 



Then, Hdt is called decoherence free if and only if pi 
undergoes p\ = — i[Hi,pi]. When we assume ^2(0) = 
P3(0) = 0, the necessary and sufficient condition for Hdf 
to be decoherence free is that the Lindblad operators and 
the Hamiltonian satisfy 



P k =a k I, Rk = 0, 



ff 2 + i$> 

k 



Qk = O, 



(5) 



where are arbitrary complex scalars [Tl|. (If pi(f)) 
and/or ^3(0) do not vanish, in addition to 10), Qk = 
O has to hold.) We then observe that the DFS always 
includes pure steady states. Actually, a pure state \tp) — 
\h) © |0) € Hdf © W^f, where \h) is an eigenvector of 
ffi and |0) is the zero vector, is a steady state of the 
dynamics because the unitary evolution pi = — i[Hi,pi] 
has a steady state p\ — \h)(h\. It is also observed that 
the matrix iff + (1/2) J^. b\.bk is now represented by 



iffi + i£ fc K| 2 iH 2 + ±J2 k a* k Q k 

[H 3 + \ J2k(QkQk + R l R k) 



2 ^k 

o 



and thus, the condition of Theorem 1 indeed holds; we 
have bk\ip) — <Xk\ip) and 



ih - 



Ctkl 



where /igl denotes the eigenvalue of ffi corresponding 
to | h) . The above discussion indicates that the condition 
for a DFS is more general than that for a pure steady 
state. 



B. The modified master equation via feedback 

The dynamical evolution of a quantum state under ho- 
modyne measurement is described by the stochastic mas- 
ter equation 0, HI El, 



dp c 
dy 



= -i[H,p, 
Tr[(c 



}dt 



\-V[c]p c dt- 
}dt + dW, 



H[c] Pc dW, 



where p c denotes the quantum state conditioned on the 
measurement data y(t), and dW is the standard Wiener 
increment with mean zero and variance dt. The super- 
operator H[c]p := cp + pc* — Tr [(c + c))p]p represents 
the stochastic jump associated with the continuous mea- 
surement. It must be noted that we consider only one 
decoherence term T>[c]p c here, which is caused by the 
interaction between the system and measurement appa- 
ratus. 

In the theory of Markovian feedback, adding a data- 
dependent feedback Hamiltonian I(t)F = (dy/dt)F, we 
can modify the properties of the decoherence as follows: 
The infinitesimal increment of the conditioned quantum 
state is given by p c {t + dt) = e I(t ^ dt!C [p c (t) + dp c (t)}, 
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where K.p = — i[F, p] . This leads to the following modified 
stochastic master equation: 

d Pc = -i[H, p c ]dt + V[c]p c dt + H[c]p c dW 

-i[F, cp c + p c c)]dt + V[F]p c dt - i[F, p c }dW. 

The controlled master equation is obtained by simply 
dropping the stochastic term in the above equation as 
follows: 



dp 
~dt 



H 



1 



-c 1 F ■ 



1 



;Fc, p 



+ V[c - iF]p. (6) 



Hence from Theorem 1 , a pure state | tf>) becomes a steady 
state of the modified master equation 10 if and only if 
the two matrices \H' + (l/2)c''c' and c' have a common 
eigenvector \<f>), where H' := H + (l/2)dF + (l/2)Fc 
and c' := c — iF. A direct calculation yields the following 
objective condition. 

Theorem 2. The modified master equation 0) has a 
pure steady state if and only if two matrices 



A = \H + \Fc + -c^c + -F 
2 2 



have a common eigenvector 
is given by p = \ip}((j)\. 



B = c-iF 
Then, the steady state 



III. EXAMPLES 
A. Purification of a single atom 

As a simple example, let us reconsider the control prob- 
lem for a single two-level atom 0. The Lindblad opera- 
tor given by 2 <J - = ^-{a x — io~ y ), where o~i ii — x, y, z) 
denote Pauli matrices, represents spontaneous emission. 
The feedback Hamiltonian and driving Hamiltonian are 
given by F = \o- y and H = aa y , respectively, where 
A £ 1 and a£l are controllable parameters. Without 
the feedback control (i.e., A = 0), only the ground state 
\g) = [0 1] T can be pure when a = 0. Let us compute 
the condition of the control parameters such that the 
modified master equation has a pure steady state. Since 
dim7Y = 2, the condition in Theorem 2 is equivalent to 
det[A, B] = 0, which leads to 



A 



Hence, a and A must satisfy the relations 

7 \F~i 

a = — sin 9 cos 9, A = — — — (l±cos0), 
4 2 v ' 

where 9 £ R is a real parameter. Then we have steady 
states 



(7) 



1^ 



cos(0/2) 
±sin(0/2) 



\4>- 



sin(0/2) 
±cos((9/2) 



where <j>± corresponds to the sign in A. It should be noted 
that Eq. Q is a necessary and sufficient condition for 
the master equation to have a pure steady state^whereas 
it was derived as only a necessary condition in 



B. Dynamical creation of a maximally entangled 
state 

The study Q successfully applied the Markovian feed- 
back method to a two-qubit system for the following 
setup. The Lindblad operator associated with the mea- 
surement is 



dy^O^ ® I + I I 




10 
10 
110 



The feedback and the driving Hamiltonian are given by 

F = \J X , H a = aJ x , (8) 

where A and a£l are the controllable parameters, and 
J x := (a x ® i" + I ® <r x )/2. Then, it was shown that the 
entanglement, which is measured by the so-called "con- 
currence" ' 01 the steady state of the modified master 
equation was larger than that without feedback. How- 
ever, the steady state was a mixed state for any of the 
parameters a and A. 

Therefore, assuming that we can make any kind 
of feedback Hamiltonian F, our aim is to character- 
ize F such that the maximally entangled state |4>) = 
(l/\/2)(|00) + 1 11)) becomes a steady state of the master 
equation. The driving Hamiltonian is given by H a in Eq. 
©. First, from the condition = (c-iF)|$) = 
in Theorem 2, the feedback Hamiltonian F = {/y} must 
satisfy 

fll = /44, J 3 4 = ~f*3 - 2^7, 

/ 24 = -A* 2 - 2^7, /i4 = -M-/n, (9) 

where we have defined p := —ik £ R. Due to the relations 
©, the condition A|$) = in Theorem 2 becomes 

V7(Aa + As) + 4 7 + M 2 /2 
ia + V7(/22 + j*2 3 ) + v^M 
ia + 77(/33 + / 2 * 3 ) + \flV 
-V7(/i2 + /i 3 )-4 7 + ^ 2 /2 

Thus, from the second and third entries, we obtain 
2ia + ^(/ 22 + fa + / 23 + / 2 *s) + 2V7M = 0, 

(/22 - / 33 ) + (f23 - / 2 * 3 ) = 0, 

which implies that a = 0, / 22 = / 33 , and / 23 £ R. Ac- 
cordingly, we have / 23 = — / 22 — p. The first and fourth 
entries in Eq. (UUJ) yield /i2+/i 3 +4 1 /7 = and p 2 = 2v, 
As a result, the feedback Hamiltonian that enables |$) 
to become a steady state is completely parametrized by 

F = (xt + x 2 )I <& I — (p + xi + x 2 )cr x ® (Tj, 
+ (x 1 — x 2 )(cTy ® cr y + a z ® a z ) 
+ («3 + n/7)U ® O - (sea + 3x/7)K <8 /) 
+ x 4 (crj / ®a z -a z ® a y ) 

- V7(o"x <8> &z + &z ® o"x), (11) 





" 1 ■ 







= ^ 







. 1 . 



(10) 
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where xi (i — 1, ... ,4) are real parameters. Especially 
when xi — X2 — x^ = p = and x 3 = the 
Hamiltonian (flip reduces to 



O ' 7 



I® Ox) 

f a z <g) a x ). 



Therefore, we need to add the global Hamiltonian repre- 
sented by the second term to the local Hamiltonian F a . 



IV. CONCLUDING REMARKS 

Although we could have found feedback Hamiltonians 
in the two typical examples such that the modified master 
equations have a pure steady state, the applicable cases 
are limited. In fact, we cannot find such a convenient 
Hamiltonian in the following two cases even under some 
ideal assumptions of the feedback scheme (e.g., negligible 
time delay of the feedback). The first is the case where 
there are uncontrollable decoherence effects bk and they 



have no common eigenvector. In this case, obviously, we 
cannot find a common eigenvector of bk and c — \F for 
any Hamiltonian F. The second case is as follows: If the 
measurement efficiency, denoted by 77, is less than 1, the 
modified master equation takes the form 



dp 
It 



-V[c-iF]p 



V 



1 



>1 



n 



(12) 



From Theorem 1, the dynamics H12|) has a pure steady 
state \(j>) if and only if \(f>) is an eigenvector of F; this 
directly indicates that \4>) is a common eigenvector of 
c and iH + (1/2)^0. Hence, if c and iH + (l/2)c t c do 
not share a common eigenvector, we are unable to achieve 
our objective. In other words, the feedback does not have 
any ability to produce a steady state. This leads us to 
conclude that we must exploit basic designing methods 
of a feedback Hamiltonian in order to obtain the purest 
possible mixed steady state. 



[1] G. Lindblad, Commun. Math. Phys. 48, 119 (1976). 

[2] D. A. L idar, A. Shabani, and R. Alicki, 

quant-ph/0411119 (2005). 
[3] M. A. Nielsen and I. L. Chuang, Quantum Computation 

and Quantum Information, Cambridge University Press, 

Cambridge (2000). 
[4] V. P. Belavkin, J. Multivariate Anal., 42, 171 (1992). 
[5] V. P. Belavkin, Theory Probab. Appl., 38, 573 (1993) 
[6] H. M. Wiseman, Phys. Rev. A 49, 2133 (1994). 
[7] J. Wang and H. M. Wiseman, Phys. Rev. A 64, 063810 

(2001). 

[8] J. Wang, H. M. Wiseman, and G J. Milburn, Phys. Rev. 



A 71, 042309 (2005). 
[9] D. A. Lidar, I. L. Chuang, and K. B. Whaley, Phys. Rev. 

Lett. 81, 2594 (1998). 
[10] D. A. Lidar, D. Bacon, and K. B. Whaley, Phys. Rev. 

Lett. 82, 4556 (1999). 
[11] A. Shabani and D. A. Lidar, Phys. Rev. A 72, 042303 

(2005). 

[12] W. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 
[13] L. K. Thomsen, S. Mancini, and H. M. Wiseman, J. Phys. 
B, 35, 4937 (2002) 



